New conformally flat initial data for spinning black holes 
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We obtain an explicit solution of the momentum constraint for conformally flat, maximal slicing, 
initial data which gives an alternative to the purely longitudinal extrinsic curvature of Bowen and 
York. The new solution is related, in a precise form, with the extrinsic curvature of a Kerr slice. We 
study these new initial data representing spinning black holes by numerically solving the Hamiltonian 
constraint. They have the following features: i) Contain less radiation, for all allowed values of 
the rotation parameter, than the corresponding single spinning Bowen- York black hole, ii) The 
maximum rotation parameter J/m^ reached by this solution is higher than that of the purely 
longitudinal solution allowing thus to describe holes closer to a maximally rotating Kerr one. We 
discuss the physical interpretation of these properties and their relation with the weak cosmic 
censorship conjecture. Finally, we generalize the data for multiple black holes using the "puncture" 
and isometric formulations. 
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I. INTRODUCTION 

Black holes are expected to be common objects in the 
universe. At the classical level, they are described by the 
Einstein field equations. The study of black holes by an 
initial value formulation of Einstein's equations is a diffi- 
cult problem, both from the analytic and numerical point 
of view. One of the most important open question regard- 
ing black holes is the weak cosmic censorship conjecture 
(cf. 0] and also the recent review |^): generic singulari- 
ties of gravitational collapse are contained in black holes. 
The physical relevance of the concept of black holes de- 
pends on the validity of this conjecture. However, a gen- 
eral proof of the cosmic censorship conjecture lies outside 
the scope of present mathematical techniques. Thus, it 
is only possible to prove it in some very restrictive cases 
such spherical symmetry or to find indirect evidence of 
its validity. In the obtaining of these indirect evidences, 
a key role is played by the following three consequences 
of the weak cosmic censorship and the theory of black 
holes i [|: 

(i) Every apparent horizon must be entirely contained 
within the black hole event horizon. 

(ii) If matter satisfies the null energy condition (i.e. if 
Tabk^k'' > for all null k"), then the area of the 
event horizon of a black hole cannot decrease in 
time. 

(iii) All black holes eventually settle down to a final 
Kerr black hole. 

From (i)-(iii) we can deduce the Penrose inequality |^ 

A < leW, (1) 

where A is the area of the apparent horizon and m is 
the total ADM mass of the space-time. Remarkably, the 



inequality (|l]) involves only quantities which can be com- 
puted directly form the initial data. After considerable 
effort, the Penrose inequality has been proved for time 
symmetric initial data |^ Q , providing an important sup- 
port for the validity of (i)-(iii). The Penrose inequality 
can be strengthened if we assume axial symmetry and 
take into account angular momentum. Angular momen- 
tum is a conserved quantity in axially symmetric space- 
times, since it can be defined by a Komar integral (cf. ^ 
and also [^). Using this fact and (i)-(iii) we deduce the 
following inequality ||l^ 



SA < 1, 



87r(m2 + Vm* - J2) ' 



(2) 



where J is the total angular momentum of the space- 
time. The inequality (^ must hold for every axially sym- 
metric, non singular, asymptotically flat initial data. The 
equality in (|^) must be achieved if and only if the data 
are slices of Kerr. Note that (0) implies 



J 



(3) 



We can also obtain upper bounds for the total amount 
of gravitational radiation emitted by the system. The 
total energy radiated is given hy E = m — mf where m f 



we have 



E 



< 1- v^. 



(4) 



In contrast with (||) and (||), the inequality (jj) involves 
the complete evolution of the initial data. 

No analytic proof is available so far for (||). On the 
other hand if the inequality (^ fails for some initial data 
then, these data represent a counter example of weak 



2 



cosmic censorship. Beside the trivial example of Kerr 
initial data, inequality (^) has been only studied in the 
spinning Bowen-York initial data (cf. [p^). 
The spinning Bowen-York data are conformally flat 
and the second fundamental form is a explicit solution 
of the momentum constraint which contain the angular 
momentum of the data as a free parameter. The mass of 
the data has to be computed numerically by solving the 
Hamiltonian constraint. Remarkably, these data reach a 
limit in and e,/ 



0.812 <eA<l, ej< 0.928. 



(5) 



eA reaches the upper limit only in the nonrotating case, 
since in this case the Bowen-York data reduce to the 
Schwarzschild data. ej cannot reach the limit case 
1 because they are not slices of Kerr for any choice 
of the free parameter J, even when J goes to infin- 
ity. It appears that the Kerr metric admits no confor- 
mally flat slices (in fact, in |15 it has been shown that 
there does not exist axisymmetric, conformally flat folia- 
tions of the Kerr spacetime that smoothly reduce, in the 
Schwarzschild limit, to slices of constant Schwarzschild 
time) . It is rather easy to construct conformally flat data 
with smaller ej: take the conformal second fundamental 
form of Bowen-York K'^y ^^nd add a solution of the mo- 
mentum constraint K'^ such that it has no angular mo- 
mentum and such that the square of the sum K'j^y + 
is bigger than the square of i^^y (these solutions can be 
exphcitly constructed using Theorem 14 of [^). Solve 
the Hamiltonian constraint with respect to the new sec- 



ond fundamental form K^^y 



K"^. Then the new data 



will have bigger mass and equal angular momentum than 
the Bowen-York one. 

It can be proved that it is not possible to reverse this 
argument to produce a data with bigger ej. Because of 
this, one is tempted to believe that (|^) is the upper limit 
to all conformally flat initial data. We will see that this 
is not the case. 

In order to test inequalities (||)-(|1) in a sharper way 
than with the Bowen-York data we need to construct 
data with higher ca and ej. To do this, it is natural to 
consider perturbations of the Kerr initial data. However 
there are many ways to perturb the Kerr initial data, and 
for each case we have to deal with the solvability of a non- 
linear elliptic equation in order to get a solution of the 
constraints. We have found a remarkable simple way of 
solving this problem. In this article we describe the con- 
struction of a conformally flat initial data in which the 
second fundamental form is related to the Kerr second 
fundamental form in a simple way. The second funda- 
mental form will be an explicit solution of the momentum 
constraint with the following property: It is a conformal 
rescaling of the Kerr second fundamental form. These 
new data can be interpreted as a conformally flat defor- 
mation of the Kerr initial data. We find that for the new 
data 



We also find that the total energy radiated is smaller 
than the Bowen-York one and satisfies the inequality (jj). 

The initial data presented here is also relevant for the 
binary problem. In order to improve the reliability of the 
results found in jl^ it is necessary to explore other family 
of initial data to know whether these results depend or 
not on the specific data used; that is, whether there exist 
physical properties of the wave form emitted by a binary 
system that are invariant under small changes on the 
data. It will be possible to measure only this kind of 
properties. Recently, other families of initial data has 
been suggested (see for instance Ref. |18 [1^ and |2^ 
PH for a Kerr-Schild type). The data constructed here 
is as simple as the standard Bowen-York, they will be a 
good candidate for future test and comparison with other 
initial data. 

In Sec. |l| we construct a explicit solutions of the mo- 
mentum constraint for conformally fiat metrics such that 
they are conformal to the Kerr second fundamental form. 
These solution are constructed in an invariant way, they 
will depend on a scalar function uj that can be explicitly 
calculated for the Kerr initial data. In Sec. [II we de- 
scribe the numerical computations of sequences of initial 
data and study the ej and ca for increasing values of J . 
We then numerically evolve these initial data to compute 
the total gravitational energy radiated to infinity in or- 
der to compare it with the Bowen-York solution (and the 
vanishing Kerr data). Finally in Sec. |l^ we discuss the 
extension of our solution to the binary (an multi) black 
hole case. 



II. CONSTRAINT EQUATIONS WITH AXIAL 
SYMMETRY AND THE NEW DATA 

The standard conformal method for solving the con- 
straints equations for maximal initial data is the follow- 
ing (cf. and the reference given there). We give 
a conformal metric hab and a symmetric trace free tensor 
if"'' such that 



ab 



0, 



(7) 



where Da is the covariant derivative with respect to hab- 
Then, we solve the following equation for the conformal 
factor (fi 



K'^'^Kab 



(8) 



0.813 <eA<l, ej< 0.932. 



(6) 



where Lh = D'^Da — R/S, R is the Ricci scalar of the met- 
ric hab and the indexes are moved with hab- The physical 
fields defined by Kb = 'p'^hab and = ip-^^K''^ wfll 
satisfy the vacuum constraint equations. We need to pre- 
scribe appropriate boundary condition to equations (|^) 
and (^), we will come back to this point later on. 

Remarkable simplifications on (^ and (|^) occur when 
hab has a Killing vector ry". We will assume that 77" is 
hypersurface orthogonal and we define 77 by 77 = ri°"rf'hab- 
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We analyze first the momentum constraint (Q) (of. pi] , 
p5| and Consider the foUowing vector field S*" 



ga ^ _,abc jj^^^ £^ ^ 0, 



(9) 



where is the Lie derivative with respect 77" and eabc 
is the volume element of hab- It follows that S*" satisfies 



0, 5"?7a = 0, DaS"" - 0. 



(10) 



Using the Killing equation Di^a'Hb) — 0, the fact that 
■q"^ is hypersurface orthogonal, (i.e.; it satisfies Da^fj — 
— r/[al?ti] In?]) and equations (!□) we conclude that the 
tensor 



our case, both condition are trivially satisfied since q will 
be chosen to be zero. For the function cj we need that 
it cancel out the singular denominator sin^ in ( [1^ ) , in 
order to obtain solutions which are smooth in 0. The 
function q determine the conformal metric, for q = we 
have that the data is conformally fiat. The function to 
determines the extrinsic curvature of the data. 

To solve ( 14|) we use the following boundary condition 
(cf. H, H, 3q ) 



lim (y5 = 1, lim r(y5 = — ^, 

r — >oc ' 1 — >0 2 



(16) 



where mp is a positive constant called bare mass. The 
total angular momentum of the data is given by 



K 



ab 



(11) 



J 



1 

'8^ 



San"" dS, 



(17) 



is trace free and satisfies 
written in terms of to 



The square of K"^'' can be 



K'^^Kab = 2- 



DcuD^u 



(12) 



Two facts are important. First, the function lo is ar- 
bitrary. In particular it does not depend on the metric 
hab- Second, the extrinsic curvature of the Kerr initial 
data (in the Boyer-Lindquist coordinates) has the form 
(prj), and then it has a corresponding function ljk- 

We can describe now the new data. Let h be the flat 
metric. Take the function ojk- Define if"* by (^ and 
(prj). Solve (^ for the conformal factor with the appro- 
priate boundary condition. We will obtain a conformally 
fiat data with a extrinsic curvature that 'resemble' the 
Kerr extrinsic curvature. In other words, we use the func- 
tion w to construct out of the Kerr extrinsic curvature a 
explicit solution of the momentum constraint (^ for flat 
metric. 

In order to write the equations and the boundary 
conditions explicitly we introduce spherical coordinates 
(r, 0, (/)), and write the metric in the form p6| 



where S is any closed two-surface which enclose the ori- 
gin. This equation can be written in a remarkable simple 
form in terms of lo 

J = i(w(r,6l = tt) - w(r,6l = 0)). (18) 

Let us mention some examples. The spinning Bowen- 
York@ initial data is obtained as a solution of equation 
(|l|) with q = 0, mo an arbitrary constant, and lo given 

by 



iOBY — J{C0S'^ 9 — 3 COS I 



(19) 



The Kerr initial data is obtained as a solution of (|l4|), 
where the function q is given by 



-29J, 



' BL 



27nKa- tbl sin 



where 



^'bl + «^ cos^ 9, 



(20) 



(21) 



-2q 



r^d9^) 



2 • 2 

r sm 



Then, the Hamiltonian constraint (^ reads 
'^4'^ 4r6 sin^ 9ip'^ ' 



(13) 



(14) 



where A is the fiat Laplacian in the spherical coordinates 
(r, 9, (j)) and 5 is the two-dimensional Laplacian 



tbl = r + rriK 



2 2 
4r ' 



(22) 



and rriK and a are the Kerr parameters, and tbl is the 
usual Brill-Lindquist radial coordinate. The function uj 
is given by 



rriKa^ sin'' 9 cos ( 



(23) 



1 



1 



Sq = -dr{rdrq) + -^d^q. 



(15) 



The Killing vector is given by 77" = {d/d(f>)°', rj — sin^ 9 
and £riLO = d^Lo = 0. Note that q and lo are almost free 
function. For q we impose that is regular at the axis. 
Also, in order to have solutions for equation (14), we 
need to impose some global condition on q (cf. [27|). In 



where, J — niKCi- And the parameter toq in Eq- ( P^ ) is 
given by 



Too ~ \/ — ■ 



(24) 



We see that the second fundamental form we are con- 
sidering can be interpreted as a correction 0{J^) (and 
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higher) to that of Bowen and York as is directly reflected 
in the form of ujx- 

The new data is g = 0, and lo and toq given by Eqs. 
(^) and ( p^ ) respectively. The explicit components of 
the conformal second fundamental form are given by 



sin^ 6*, 



(25) 



-2a^mKrBL 



4r^ 



cos f sm 



Note that mo is chosen such that we recover the Kerr 
initial data if we insert in equation ( p^ ) the corresponding 
function qK- 

Since g = 0, in this case the boundary condition jl^ ) 
can be achieved by the ansatz 0=1 + mo / (2r) + u, where 
the function u is bounded at the origin and satisfies the 
regular equation in 'E? 



Au = 



(26) 



with the boundary condition lim^— >oo u — 0. This allows 
to use the 'puncture' treatment in numerical implementa- 
tions of these initial data [^ . Existence and uniqueness 
of positive solution of the elliptic equation (|6|) can be 
proved by standard methods (see for example and 
| p6[ and reference therein). 

Finally, we want to prove that these data have an isom- 
etry defined by an inversion map through a sphere of ra- 
dius TOo/2- Inversion transformation in the context of 
black holes initial data have been studied in jl^. The 
inversion map / is defined in spherical coordinates by 



(27) 



We want to prove the following assertion: if the functions 
uj and q satisfy 



I o q = q^ I o LO ~ UJ, 



(28) 



(that is, they are invariant under /); then the physical 
fields satisfy 



{rh)ab = Kb, {I*K)ab = -Kab 



(29) 



In the examples presented here, one easily check that lu 
and q satisfy Eq. (|2^), since their radial dependence is 
given in terms oi tbl defined by Eq. (p2[), which satisfy 
I ° rsL = rsL- 

In order to prove the assertion, we first calculate the 
transformation of the conformal fields under /, from Eqs. 
( prj ) and (|l^) we obtain 

(30) 
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FIG. 1: Rotation parameter ej as a function of the angular 
momentum J for the new data, Kerr, and Bowen-York holes 
along the curve, keeping fixed the free parameter mo = 2. 
The curve for Kerr can be obtained analytically by equation 



Then, using Eq. (||) we obtain 

I o in ^ { — ] 



\2r J 



Lp. 



(31) 



From Eqs. (gOj) and (|3|) follow Eqs. (|29|). This property 
of the data is used in the numerical calculations, as we 
will see in the following section. 



III. NUMERICAL RESULTS 

The 2D fully nonlinear evolutions have been performed 
with a code, Magor, designed to evolve axisymmetric, 
rotating, highly distorted black holes, as described in 
Ref. [g^; |3^. This nonlinear code solves the complete 
set of Einstein equations, in axisymmetry, with maximal 
slicing, for a rotating black hole. The code is written 
in a spherical-polar coordinate system, with a rescaled 
logarithmic radial coordinate that vanishes on the black 
hole throat. An isometry operator is used to provide 
boundary conditions on the throat of the black hole. The 
lapse is chosen to be antisymmetric across the throat. 
The shift vector components are employed to keep off 
diagonal components of the metric zero, except for the 
gs^. The initial data described in Sec. |ll| above are pro- 
vided through a fully nonlinear, numerical solution to the 
Hamiltonian constraint. 

In Fig. |l| we plot the curves of constant bare mass 
mo = 2 for initial data corresponding to rotating holes of 
Kerr, Bowen-York and the new ones introduced in this 
paper. Interestingly enough it was shown in Refs. p^ , 
|l3| that the Bowen-York hole reach a maximum rotation 
parameter of cj w 0.928 when J goes to infinity. We 
have reproduced this curve as a function of J reaching 
values of J w 10, 000. For Kerr holes, along the curve 
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FIG. 2: The horizon area parameter as a function of the 
angular momentum J for the new data, Kerr, and Bowen- 
York holes along the curve, keeping fixed the free parameter 
mo = 2. For the Kerr initial data the exact value is = 1. 

Too = constant we have equation (|2^); the extreme Kerr 
a/niK — 1 correspond to the limit J — > oo. We find for 
the data presented in this paper that the maximum lies 
near ej = 0.932, which is higher than the Bowen-York 
maximally rotating hole. This allows to study black hole 
evolutions from values closer to the maximally rotating 
Kerr ones. 

The code is able to evolve such data sets for time 
scales of roughly t < IOOto, and study physics such as 
location and evolution of apparent horizons and gravita- 
tional wave emission. We have used typical grid sizes of 
300 radial by 39 angular zones and extracted waves at 
the radial location robs = 15to. The results are displayed 
in Fig. ^ and clearly show that the new data has less 
radiation content than the spinning Bowen-York holes. 
It produces roughly ten percent less the total radiated 
energy, E. Also this plot shows that inequality (^) is 
satisfied, the upper bound given by (Q) is in this case 
w 0.031 while the maximum of the total energy radiated 
is w 0.0015. 

For completeness we have plotted the evolution of Kerr 
initial data, for which the outgoing radiation should 
strictly vanish, as measure of the numerical error of 
our evolutions; and for comparison with future work we 
present the numerical parameters of our initial data fam- 
ily in Table | 

IV. DISCUSSION 

On the light of the two aspects studied in this paper we 
observe that the new data improves on the Bowen-York 
one leading to less spurious radiation and allowing a rep- 
resentation of higher rotating black holes while keeping 
the simplicity of the solutions, namely the explicit ana- 
lytic form of the extrinsic curvature and conformal flat- 
ness of the three-geometry. This proves that even within 




FIG. 3: Radiation content of a single rotating black hole 
given by the new, Kerr and Bowen-York initial data. The 
upper bound given by (0) is in this case « 0.031. 



TABLE I: Initial spin J, ADM mass m, apparent horizon mass 
niAH = \/ A/IGtt, and energy radiated of the new family of 
initial data 



J/(mo/2)' m/(mo/2) niAH /{mo/2) E/m, 



1 


2.0478 


2.047 


2.75e-6 


2 


2.1694 


2.169 


2.99e-5 


5 


2.674 


2.668 


0.00028 


10 


3.475 


3.452 


0.00066 


20 


4.743 


4.686 


0.00100 


50 


7.375 


7.242 


0.00131 


100 


10.386 


10.165 


0.00140 


500 


23.169 


22.580 


0.00162 


1000 


32.760 


31.899 


0.00163 


5000 


73.245 


71.248 


0.00191 


10000 


103.583 


100.74 


0.00187 



the conformally flat ansatz one can look for astrophysi- 
cally more realistic initial data. Also the new data satis- 
fies inequalities (||), (||) and (||), supporting the vahdity 
of weak cosmic censorship. 

We want to discuss now the generalization of these data 
for multiple black holes. In general, if the conformal met- 
ric hab admits only one axial Killing vector 77" the only 
freedom left is the choice of the origin in the z = cos6/r 
coordinate. By superposing different solutions of the mo- 
mentum constraint of the form (|l^) such that they are 
singular at different points in the axis we will obtain a 
multiple black hole solution for a general axially symmet- 
ric metric. The spin of all of the black holes will point 
in the direction of the axis. However, when hat is cho- 
sen to be the flat metric, then we have that a rotation 
about any axis is a Killing vector. Hence it is straight- 
forward to generalize the data presented here to include 
multiple black holes in arbitrary location and with spin 
pointing in arbitrary direction. For completeness we will 
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where Cartesian coordinates, x"^ and z° are arbi- 

trary constant vectors (we chose z° to be a unit vector) 
and — (x'^ — x'^){xc — Xc). The vector x'^ represent the 
new origin and the new axis. The new coordinate 9' 
with respect to the axis is given by 



cosf 



Zcix" - x"). 



(33) 




FIG. 4: Geometry in the conformal (Cartesian) space of the 
location and orientation of a generic nth spinning black hole. 



Fig. U displays explicitly those vectors. 

Then, the desired expression for K"-^ is given again by 
Eqs. (prj), (||) and ( p3| ) but we use in equations ( |ll| ) and 
(^) the expression for 77" and 77 given by (^^, and we 
replace in (|2^) 9 by 9' given by (33). 



write explicitly the expression for the general solution of 
the momentum constraint, for flat metrics, with arbitrary 
origin and with spin pointing in arbitrary direction. The 
general Killing vector for the flat metric can be written 

as 

r;"^ = e''^^z,(x,-xe), 11 ^ - {z,{x' - x'))\ (32) 
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